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Abstract 

We show that a model, recently introduced for quantum nondemolition measurements 
of a quantum observable, can be adapted to obtain a measurement scheme which is able to 
slow down the destruction of macroscopic coherence due to the measurement apparatus. 


1 Introduction 

One of the most important limitations in the observation of quantum coherence at macroscopic 
level is the possibility of generating at least to macroscopic quantum states which show the quan- 
tum coheh rence. Since the seminal work of Yurke and Stoler [1] it becomes clear that a Kerr 
medium could be used to generate such states at optical level. They showed, indeed, that the 
unitary evolution of an initial coherent state, interacting with a Kerr medium with a well defined 
length, will produce a superposition of coherent states. For instance an initial states |a) will 
generate the superposition 

l (e -^ |a) + e W Va)) (1) 

after an interaction time t 0 = tt/(2D) where D is the strenght of Kerr nonlinearity. At well defined 
shorter times three or more coherent states could also be generated [1]. This, of course, requires 
the precise knowledge of the length of the medium (or interaction time). It is also well known, 
and was shown in great details by Daniel and Milburn [2], that as soon as one takes into account 
the loss in the Kerr medium the generation of those states is suddenly inhibited. Thus, the best 
should be to have a Kerr medium with high nonlinearity to loss ratio. Recently [3], quadrature 
squeezed light was observed in semiconductors at frequencies less than half of band gap, where 
large ratios of nonlinearity to loss can be obtained [4]. Then, semiconductors could be the best 
media to generate the superposition of states because of the large ratio of the nonlinear phase 
shift to the optical losses which in the reported experiment [3] was extimated greater than 100. 
Furthermore, it has been recently shown [5] that a quasi-superposition of macroscopic states, with 
interference fringes still present, could be generated in a Kerr medium with the above ratio of 10, 
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when one uses a squeezed bath to model the loss. In this context it was also shown that a squeezed 
bath could be realized by a suitable feedback [6]. Moreover, it was also shown [7] that by using 
a time modulation of the Kerr nonlinearity one could obtain the coherent superposition without 
the precise knowledge of the length of the medium (or interaction time) by only adjusting the 
phase of the time modulation. However, even though we could assume that such a macroscopic 
superposition (or quasi-superposition) has been generated, one should have some experimental 
apparatus suitable to observe the interference pattern. Yurke and Stoler [1] pointed out that 
any unavoidable dissipation, introduced by the measurement process, will suddenly destroy the 
interference fringes which are the signature of the coherent superposition. Kennedy and Walls [8], 
following a suggestion of Mecozzi and Tombesi [9], showed that a phase-sensitive experimental 
apparatus, like the one modeled by a squeezed bath, might preserve the macroscopic coherence. 
In the present paper we will show that such an experimental device could be phisically realized by 
using an appropriate quantum nondemolition (QND) model, introduced by Alsing, Milburn and 
Walls [10], when one takes into account the detunings of the coupled modes with respect to the 
cavity characteristic frequencies. 
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2 The Model 


We consider a cavity supporting two different modes, with annihilation operators a and b. The 
two modes are coupled by a nonlinear crystal, so that (in the interaction picture) 


Hint = h X X e Y v , 


( 2 ) 


where X £ — (ae* + a f e /2 and Y v = (be' v + tfe /2. This interaction could be achieved 

by, for example, a crystal with a x ^ nonlinearity in which two processes driven by classical 
fields, amplification at the frequency u, = u a + Ub, and frequency conversion at the frequency 
u>d = u a — Ub, have equal strengths [10]. Because of the QND condition, when the “meter” mode 
b is heavily damped at rate kb, one can monitor the quadrature of the signal mode a just by 
performing a homodyne measurement of a quadrature Y f of the mode b. In fact, when kb » k a 
(damping rate of the a mode) the homodyne photocurrent I(t) can be directly expressed in terms 
of the “instantaneous” mean value (X^(t)) c , conditioned on the result of the measurement [11, 12], 
as 


I(t) = 7]X 


2 sin (5 -ip) {Xf.(t)) c + 



( 3 ) 


where rj is the efficiency of the homodyne detection and £(f) is a Gaussian white noise with 

The QND-mediated feedback model of [6, 1 1] is obtained by taking part of the output homo- 
dyne photocurrent I(t) and feeding it back to the cavity so to add a driving term H/b(t ) = hgI(t)X 0 
to the a mode Hamiltonian. The constant g represents the gain of the feedback process and 
X 0 = (ae* 0 + a^e~' e 'j /2. If one adiabatically eliminates the meter mode b and applies the Marko- 
vian feedback theory recently developed by Wiseman and Milburn [13], the dynamics of the a 
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mode can be exactly determined, and in [11] we have shown that in the unstable regime the 
decoherence time of an optical Schrodinger cat can be appreciably increased, so to facilitate its 
detection. 

In the present paper we reconsider this model and we eliminate the electro-optical feedback 
loop. We simply detune the two modes in the cavity, so that their uncoupled evolution is no more 
driven by the standard vacuum bath term alone, but by 


C a p = k a (2 apa S — a)ap — pa^aj — i [^a^p] 


( 4 ) 


and an analogous expression holds for the b mode. The effect of the two nonzero detunings S a 
and 6 b can be intuitively described in terms of an “internal feedback” mechanism, because the 
detunings mix the two quadratures X c and Y v with their respective 7t/2 out of phase quadrature, 
so that any variation of X; is “fed back” to the X ( dynamics itself by the joint action of the 
detunings and the nonlinear coupling. Provided that the adiabatic condition k b k a is satisfied, 
the homodyne measurement of the quadrature Y s allows monitoring the a mode quadrature X ( 
also in the presence of nonzero detunings. In fact, when 6 b ^ 0, Eq. (3) generalizes to 


I(t) = TJX 


2k > ' ( f\ ^ 

WTs! sm(6 - v) 


2k b 6 b _ 
k b + <5 b 


cos (6 — <p) 




( 5 ) 


so that from the homodyne photocurrent it is still possible to reconstruct the marginal probability 
distribution of the quadrature X c , which is the quantity usually considered for revealing the 
interference fringes associated to an optical Schrodinger cat. We have therefore the model defined 
by the following master equation for the density matrix D of the two modes 


D = C a D + C b D-UH int ,D] , (6) 

h 

where the superoperator C x (i = a, b) is given by (4). We shall now see that all the interesting 
results obtained for the feedback model of [11] (the preservation of macroscopic quantum coherence 
in particular) can also be obtained with this simpler model. 

Eq. (6) can be exactly solved, because the Wigner function of the two modes evolves according 
to the Fokker-Planck equation for a four-dimensional Ornstein-Uhlenbeck process [14]. Anyway, 
the analytical expressions in the general case are very cumbersome and therefore we shall ex- 
plicitely discuss only the adiabatic limit k b » k a , where the meter mode b can be adiabatically 
eliminated, and which, as we have seen above, is the most interesting case for our purposes. After 
the adiabatic elimination of the b mode, one gets the following master equation for the a mode 
reduced density matrix p 

i> = c«p- 5 [^[x ( ,p]]+«f’ [*«,{*(,*>}] . m 


where T = + il), F = x 2 &/4(kf + ©. 
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3 Macroscopic Coherence 

We will now focus on the detection of optical Schrodinger cats rather than on their generation, and 
therefore we shall assume that at t = 0 a superposition of coherent states of the a mode has been 
already prepared, i.e. , we consider an initial condition p(0) = H Q ,p N a< p\a) (J3\. The exact time evo- 
lution from this initial state can be obtained with the same method of [1 1] and it is better expressed 
in terms of the normally ordered characteristic function x(A, A*; t) = TV |p(t) exp(Aa^) exp(— A*a)| 

x(A, A*; t) = J2 N a ,M<x) exp {S*(t)A - A (t)A* (8) 

a,/3 

-H<)|A | 2 + ^a - 2 + ^a 2 }, 

where 


A(t) = 


a 2A — iF + 2 i6 a ^ iFe 2< * 


2A 


4A 


+ 


a 2A + iF — 2 i6 a ^ „ iFe 2 * 


2A 


4A 


,-(k a +A)t 




(9) 


B*(t) = 

+ 


[3* 2A + iF- 2 iS a iFe 2 * 

ha 


2A 


4A 


(3* 2A - iF + 2i8 a iFe 2 * 
a 


2A 


4A 


(fco+A)t 


a -(k a -A)t 


(10) 


F /T<5 a 2 F\ f l - e~ 2{ka+A)t \ 

~ 16 V A 2 A J V 2 (Jt a + A) J 
r 6 a { 28 a - F) 


8A 2 


1 — e -2 *^^ 

2 ka ) 


F /T<5 a 2FWl 
+ 16 + AT/ V 2 (k a - A) / 


( 11 ) 


_ F 2 e~ 2ii (Y6 a 2 F\ l - e ~ 2 ^ +A)t 

" 16 (2<5 a -F + 2iA) V A 7 ~~ AT J 2 (fc a + A) 
FTe- 2ti 6 a 1 - e~ 2kat 
8A 2 W a 

F 2 e~ 2ii (T8 a 2F\ 1 - e - 2( *“~ A)£ 

16 (2<5 a - F - 2iA) \ A 2 + "A J 2 (Jb a - A) 

A = \J& a F — <5 2 . 


(12) 


(13) 
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We see that the system is stable and reaches a steady state if and only if 

A < A ; tt i.e. x 2 6<A < 4 (k 2 b + <5 b 2 ) (fc 2 + 6 2 ) . (14) 

In the stable case, the stationary state is described by a Gaussian density operator of the form 

p st = Z _1 exp na^a - ya t2 - ^-a 2 J , (15) 

where Z is a normalization constant and the equilibrium parameters m and n can be written as 


n 


I'oo + 1/2 


X log< 


y ^ oT +^ y ^) — l /^ ool 2 

sj (i'oo + 1/2) — l/iool 2 + 1/2 


(16) 


VteWco + 1) “ l/*oo| 2 


m = 


/x 0 


I'oo + 1/2 


n , 


(17) 


where the asymptotic values J/o© and /-too are easily obtained from (11) and (12). An interesting 
aspect of this stationary state is that it can show arbitrary quadrature squeezing. For example, 
the stationary variance of the quadrature X € is given by 


<*l> = 5 


1 8g (rgq + 2 Fkg) 

2 + 8fc 0 (A: 2 - A 2 ) 


(18) 


and one has squeezing when S a Sb < 0 and kb/k a < |<5&/<5 a |. H Is easily seen that when 8 a 0 
no squeezing is possible, while for 8 a / 0 but 8b — 0 extra noise is added to the system. The 
possibility to obtain squeezing with this model is thus only due to the existence of detunings, 
which give a sort of implicit feedback. 


4 Interference Fringes 

Let us now focus on the detection of the interference fringes associated to a linear superposition 
of coherent states. These fringes can generally be seen from the marginal probability distribution 
of the quadrature X t , P{x t ) = (x c |/>(t)|x € ), where |x e ) is the eigenstate of X ( with eigenvalue 
x e . As we have seen above, this probability distribution can be reconstructed from the homodyne 
measurement of the meter mode b and its general expression can be easily obtained from the 
characteristic function (8) [8, 11] 


Pfaft) = r 

a,0 y‘ 


(M 


7T<7; 


■(t) 


: exp 


(- 


(s { - 




(19) 
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where 


^ + v(t) + He {/^(i)e 2 ^} , 

**«) = 


( 20 ) 

( 21 ) 


As a special case we consider the initial superposition trated by Yurke and Stoler [1], produced 
by the unitary evolution of a coherent state in a Kerr medium 


P(0) (22) 

= i (e~ i1,,A \a) + e** \ - a)) (e‘" 4 <a| + e-‘"V a l) • 

With this choice (19) simplifies to 


p fa,t) = ^ [p\{x^t) +p 2 _{xz,t) 

+2p + (x i ,t)p_(x i ,t)sin[Ct(x 0 t)] |(a| - a)!^} . (23) 

The first two terms p±{x^ 9 i) describe the two Gaussian peaks corresponding to the two coherent 
states | ± a) of the initial superposition and they are explicitely given by 


where 


p±(x ( ,t) 


' exp /_(*** * ++ W)‘ 








<5 ++ (i) = Re|ae 1? r(t)| 

r(t ) = e~ kat ^cosh At — z^-sinhAtj . 


(24) 


(25) 

(26) 


The third term in (23) describes the quantum interference between the 
the function 


= ^ylm{ae^r(<)} 


two coherent states, where 


(27) 


gives the probability oscillations associated with the interference fringes and the factor |(a| — 
a )|»7(<) = exp {— 2|a| 2 r](t)} describes the suppression of quantum coherence due to dissipation. It is 
clear that this suppression is practically immediate for macroscopically distinguishable states (i.e., 
large |o|), unless T)(t) — 0. It is therefore important to analyze the behavior of this decoherence 
function Tj(t), which is equal to 


T](t) = 1 - 


Mi! 

2am 


(28) 


To be more specific, if we want to determine the conditions under which the detection of macro- 
scopic quantum coherence is facilitated, we have to compare T)(t) with the corresponding decoher- 
ence function of a standard vacuum bath, which is given by [8] 


Vvac(t) = 1 - e 2kat . 


(29) 
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This function shows that in the standard case, after a time t ^ 1/(2 k a ), it is f] vac (t ) ^ 1 and 
therefore the quantum interference is quickly washed out. On the contrary, in the present model 
it is possible that rj(t ) assumes much smaller values, so to significantly slow down the destruction 
of the interference pattern. 


5 Conclusions 

Differently from a very large part of the literature on optical Schrodinger cats, we have focused 
on their detection rather than their generation because, as realized since the paper by Yurke and 
Stoler [1], to detect a linear superposition of macroscopically distinguishable states is more difficult 
than to create it. To the best of our knowledge, only the paper by Brune et al. [15] affords a 
detailed discussion of both aspects, our opinion, Brune large number of atoms reconstruction of 
the probability distribution revealing the contrary shows how to prepare a fully optical detection 
scheme based on a very simple model, offer a promising way to both measurements and detect a 
linear superposition of coherent states. 
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